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Abstract 

The finite-time exergoeconomic optimal performance for a quantum Stirling engine is studied. The 
maximum exergoeconomic profit, the optimal thermal efficiency and the optimal power output of the 
quantum Stirling engine are obtained. © 1999 Elsevier Science Ltd. All rights reserved. 
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1. Introduction 

Stirling engine is an external combustion engine with brilliant prospects. It utilizes the resources 
of low price and easily obtainable energy such as solar, atomic, geothermal, industrial waste heat 
and marsh-gas etc. It has very important significance for environmental protection and moder¬ 
ating the tense petroleum needs in the world. With this great potential, the Stirling engine has been 
used widely in energy resource utilization and power engineering. Some authors have studied the 
performance optimization problem for the Stirling engine using the method of finite-time ther¬ 
modynamic analysis [1,2]. Significative results have been obtained by Wu [3], Hu [4] and Ladas 
and Ibrahim [5]. However, in the work mentioned above, the classical engine with pure ther¬ 
modynamic performance objective was studied while the quantum performance of engine’s 
working fluid was not considered. The dynamical semi-group approach of the quantum theory of 
open systems has been introduced into finite time thermodynamics only recently [6-12]. Based on 
the Salamom and Nitzan’s work [13] in which the economic index of profit was taken as the 
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optimization objective, Chen [14,15], Wu [16] and DeVos [17] developed the method which is 
called finite-time exergoeconomic analyses or endoreversible thermoeconomics. The exergoeco- 
nomic optimum performance of a quantum-mechanical Stirling engine is studied in the present 
paper by using the modern non-equilibrium quantum statistical theory and the methods used in 
Refs. [6-17], The maximum profit and the optimal thermal efficiency of the engine at the classical 
limit are derived. The results provide a solid theoretical basis for design and operation of Stirling 
engine. 


2. Quantum-mechanical Stirling engine 

The working fluid of the engine consists of many non-interacting harmonic oscillators or many 
non-interacting two-level systems. The former case will be carried out in this paper. The latter case 
has been analogously studied by Feldmann [11], The population of the oscillators can be obtained 
from the Bose-Einstein distribution [12,18,19], 

77 = 1/[exp(/To) - 1], (1) 

where to > 0 is the frequency of the oscillator and /] = 1 /T (T is the internal temperature of the 
working fluid). 

We assume that the engine operates between an isothermal heat source B\ at a temperature T h 
and a heat sink B 2 at T c . The endoreversible quantum-mechanical Stirling engine is made of two 
isothermal branches connected by two regenerating branches (see Fig. 1). Because of the heat 
transfer effect, the internal temperatures of the working fluid are 7) and 77 in the two isothermal 
processes, e.g. 1-2 and 3^1 in Fig. 1. The second law of thermodynamics requires 
Pc>02>P l > Pin where // = l/T t (i = h,c , 1 and 2). 

In the regenerating processes, the two pistons of the Stirling engine move simultaneously with 
the same velocity and direction. The volume of the working fluid is therefore invariant. The 
quantum working fluid flows through the regenerator and exchanges heat with the filler inside the 
regenerator. The regenerator has a high frequency and a short period. The population n of the 
harmonic oscillators increases drastically with the decrease in the reciprocal (/I) of temperature. 
Hence, the oscillator frequency in the regenerating processes is approximately constant. In Fig. 1, 
the frequencies of the two regenerating processes connecting the two isotherms are given by to = 
oj i and to = to 2 with to\ > to 2 . The cycle of the engine is illustrated in the (to, n ) plane. In the 
quantum calculation, we set h' = h/2n = 1, where h is the Planck’s constant. 


3. Cycle period 

The Hamiltonian of the harmonic oscillator system S is described in the following form [8,12] 

H s = coN = co(a + )a, (2) 

where N is the number operator, a + and a are the Bosonic creation and annihilation operators 
{N = (a + )a; [a, a + ] = 1}. 
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Fig. 1. Endoreversibe quantum Stirling cycle in ( n,co ) plane. 


In the isothermal processes, the working fluid S becomes an open system owing to the coupling 
action between S and the heat source B x (or heat sink B 2 ). Using the generalized master equation 
in an open system and in Heisenberg picture, one obtains the motion equation of the operator x in 
the system S [20] as shown in the following equation. 

dx/dt = iL f x - ^{K-) + [x, Q A ]Q y - (m 7i )'Q y [x, Q*]}, (3) 

where ( to kj ) + and (o> /k ) are constants related to the matrix element of the heat source B x (or heat 
sink B 2 ) density matrix, Q, and Q k are the operators of the working fluid system S, [,] is the 
quantum Poisson’s bracket, while 


L/ — L s + T yB h SB p B , (4) 

where Ly is the Louiville’s operator of the system S, T yB represents the trace of the operator of the 
heat reservoir (B i or B 2 ), L SB represents the coupling action between the working fluid S and the 
heat reservoir (B x or B 2 ), p B is the density matrix of the heat reservoir (B x or B 2 ). Because the heat 
reservoir (B\ or B 2 ) is very large and its internal relaxation is very strong, the disturbance of the 
working fluid on the reservoir can be neglected. Thus we obtain 


T.bLsbPh — 0 . 


(5a) 
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Using the definition of Louiville’s operator for the number operator N, one obtains 


L 5 N = [H 5 ,N] = [toN.N] =0. (5b) 

Combining Eqs. (4), (5a) and (5b) gives 

L/N = (Ly + T yB L SB p B )N = 0. (6) 

Setting x = N, Qi =a + , Q 2 = a and combining Eqs. (3) and (6) gives 

dn/dt = dN/dt = —[(coi 2 ) + + (co^) ] — [(ct> 2 i) + + (c<Ei) ]« + [(® 2 i) + + (® 2 i) ], (7a) 

where N = « is the expectation values of the number operator N, [(co 12 ) + + (^ 12 ) ] and [(co 2 i) + + 
(ca 2 i) — ] are constants. The solution of Eq. (7a) is given by 

n = n e + (n a - n e ) exp{-{[(co 12 ) + + (co 12 ) - ] - [(co 2 i) + + (co 2 i) ]}?}, (8) 


where n e = [(co 2 i) + + (co 2 i) — ]/{[(co i2 ) + + (<wi 2 ) - ] - [(® 2 i) + + (c^ 2 i) ]} is the asymptotic stationary 
value of n, n a the value of n at t = 0. Eq. (8) is similar to Eqs. (3) and (4) of Ref. [8], Using the 
same method as used in Ref. [8], one obtains [(co 2 i) + + (co 2 i)~] = a[exp (q[} -(d)] and 
[( C oi 2 ) + + (can) - ] = u{exp[(/ + q)fijCo]}. Substituting them into Eq. (7a) yields 

dn/dt = —a[exp(qpjCo)]{[exp(pj(o) — \}n — 1}, (7b) 

where a and q are constants, 7} = 1 = h, c ) are the temperatures of the heat reservoirs, n = N 

is given by Eq. (1). 

Substituting Eq. (1) into Eq. (7b), expanding the time integral and taking the second order in 
/ho and p/o at the classical limit (n » 1 and fico C 1) gives 

t\ = / [(d«/dto)/(dn/d^)]dto = A/(p l - p h ), coi < to < co 2 . (9) 


t 2 =A/(P c -P 2 ), (10) 

where t\ and t 2 are the isothermal processes times of processes 1-2 and 3-4, A = [(co 2 ) 1 — 
{of ']/« is a positive constant in matters of the quantum performance of the working fluid. 

Assuming the times of the regenerating processes 2-3 and 4-1 are t 3 and t 4 , respectively, the 
cycle period time, t, can be given by the following expression 

t = t\ + t 2 + t 2 T t 4 . (1 la) 

Let b = (t\ + t 2 )/x denote the ratio of the isothermal time to the total cycle period time, one 
obtains 


T — (h + t 2 )/b — A[(P l — p h ) + (P c — P 2 ) }/b. 


(lib) 
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4. Maximum profit 

From Eq. (2) the rate of change of energy of the working fluid is obtained 

dE/dt = d(H s )/d? = {dco/dt)n + co(dn/dt). (12) 

Comparing Eq. (12) with the time derivative of the first law of thermodynamics, we have the 
instantaneous power and the instantaneous heat flow, respectively: 

dW/dt = — (dco/dt)/?, (13) 

dQ/dt = (dco/dt)n. (14) 

Substituting Eq. (1) into Eq. (14) and expanding the heat quantity integral to first order in 
P,o)(i = /?, 1 or c, 2) gives 

r co d« = X/f5 x ,n\ < n < » 2 , (15) 


02 = X/P 2 , (16) 

where 0i and Q 2 are the absorbing heat quantity in the isothermal process 1-2 and the exothermic 
quantity in the isothermal process 3-4, respectively, and the positive constant X = ln(coi/co 2 ) is the 
natural logarithm of the frequency ratio. 

The heat exchange of processes 2-3 and 4-1 are carried out inside the engine. If the regenerator 
is free of loss, the quantity of the regeneration is independent of the energy dissipation of the 
cycle. Thus the output work and the thermal efficiency of the cycle are 


IF = 01 — 02 

(17) 

and 


n1 = (0i - 02 )/ 0 i = 1 -x, 

(18) 

where x = If/If. 



In order to obtain the optimal economic performance we must study the profit of the engine. 
The profit is defined as the difference between the revenue rate of the work output and the cost 
rate of the exergy input [13-16], The exergy input of the cycle is equal to the reversible work of the 
cycle with the same working condition. The reversible work produced by the engine with the same 
working condition is 

W r = 0j( 1 - TJT h ) - 0 2 (1 - T a /T c ) = c\Q\ - e 2 0 2 , (19) 

where T a = 1 /fi () is the ambient temperature, ei = 1 — T 0 /T h = 1 — p h /p 0 and e 2 = 1 — T a /T c = 1 — 
p c /ft 0 are the Carnot coefficient of absorbing heat and removing heat, respectively. 

We assume that the prices of exergy input and work output are cf >i and fc- The profit of the 
engine is 
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n = (02^ — 4>iW r )/x. (20a) 

Substituting Eqs. (lib), (15), (16) and (19) into Eq. (20a) yields 

n = {fbfi 2 /A)[( 1 — re\)/x — (1 - re 2 )][0 2 /(0 2 x - /?,,) + p 2 /(P c ~ 0i)Y\ ( 20b ) 

where r = <f> x /<f> 2 is the price ratio. In order to assure the engine be profitable, the price ratio r 
should satisfy the relation of 0 < r < 1. 

Using Eqs. (20a) and (20b) and the theory of extreme value of a binary function, the maximum 
profit is obtained as follows 

n m = (Xb<f> 2 /A)( 1 - re 2 )[(f - e 3 )/(l + c 3 )] 2 , (21) 

where / = [(1 — re\)/{\ — re 2 )] 112 , e 3 = (/?/,//? c ) 12 = (7^/7),) 1/2 . The conditions at the extreme 
point corresponding to Tl m are: 

/? 2o = /? c (l + e 3 )/(l+/), (22) 

x 0 =/e 3 . (23) 

The corresponding optimal working temperatures of the working fluid in the isothermal pro- 


cesses are: 


Tlo = 1/020, 

(24) 

T\ 0 = 1 /(P 2o Xo). 

(25) 


5. Discussion 

(1) Optimal thermal efficiency. The thermal efficiency of the engine that corresponds to the 
maximum profit is 


tjo=l-Xo = l- e 3 [(l - rei)/(l - re 2 )] 1/2 , (26) 

where f/ 0 is the thermal efficiency of the quantum-mechanical Stirling engine for the finite-time 
exergeoeconomic performance bound at the classical limit. 

(2) Optimal power output. The optimal power output of the engine corresponding to the 
maximum profit is 

Po = {W/x) 0 (27a) 


or 


P 0 — {Xb/A){ 1 —/e 3 )(l — e 3 //)( 1 + e 3 ) 


(27b) 
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When the power output is P„, the engine has the optimal economic benefit and its profit has the 
maximum value /7 m . 

(3) Influence of the price ratio. The price ratio r directly affects the thermodynamic and eco¬ 
nomic performances of the engine. With the typical parameters of working conditions 7), = 1000 
K, T c = 290 K, T 0 = 300 K and Xb/A = fbf 2 /A = 1, the relationships of TI m — r,P„ — r and if — r 


are shown in Figs. 2 and 3, respectively. 

When the price of the work output is far higher than that of the exergy input, r = <j) l /(f> 2 —> 0, 
Eqs. (26), (27b) and (21)) become 

ho = >1 ca, (28) 

P 0 = (lb/A)(r, CA ) 2 /(2-r, CA ) 2 , (29) 

Il m = (2bcj) 2 /A)fl CA ) 2 /(2 - rj CA ) 2 = cj) 2 P 0 , (30) 


where r\ CA = 1 — e?, = 1 — ([),,//(.) 12 = 1 - (T c /T h )' n is the C-A (Curzon-Ahlborn [21]) efficiency. 
Eq. (30) indicates that in the limit r —> 0, the maximum profit index of the engine is transformed 
into the maximum power output index and the performance bound is the C-A efficiency Oca ■ 



Fig. 2. Relationship among optimal power output (P a ), maximum profit (77 m ) and price ratio (r). 
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Fig. 3. Relationship between optimal efficiency ()/„) and price ratio (r). 


When r = 4>i/4 > 2 0, Eqs. (26), (27b) and (21) become 

i 0 = ic, ( 31 ) 

Po = o, (32) 

n m = 0, (33) 

where rj c = 1 — fi h /fl c = 1 — T c /T h is the Carnot efficiency corresponding to a classical reversible 
Stirling cycle. 

In general case of 0 < r < 1, the optimal efficiency of the engine satisfies 

Ica^Io < 1c- (34) 


(4) Influence of the quantum performance of the working fluid. In Eqs. (21), (26) and (27b) the 
constants A and X are functions of the oscillator’s frequency coi and co 2 . They depend on the 
quantum performance of the working fluid. Therefore, under given working condition, choosing a 
suitable quantum working fluid to enlarge the constant XIA is very helpful for improving the 
thermodynamic and economic performance of the engine. 
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6. Conclusion 

Economics is a major index for a real thermal power cycle. Combining finite-time thermody¬ 
namics and non-equilibrium quantum statistical mechanics with exergoeconomics, this paper 
derives the maximum profit and its corresponding performance bound (rj 0 ) of an endoreversible 
quantum Stirling engine. The result of the paper provide a profit bound for designing a real 
Stirling engine whose working fluid is a quantum fluid, such as He, Ne, D 2 etc. In the work of 
Feldman et al. [11], the cycle of operation is analogous to a four-stroke Otto cycle and the work 
output was taken as optimization objective. The major advances of this paper over Refs. [6-12] 
are: (1) the modeling of a quantum Stirling engine with working medium consisting of non-in¬ 
teracting harmonic oscillators, and (2) the introduction of a profit objective. A future step of this 
analysis is to extend the working fluid model from the harmonic oscillators or two-level system to 
a three-level amplifier as done by Geva and Kasloff [9,10]. 
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